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The basic thermodynamic parameters of the superconducting state in YNi 2 B 2 C were calculated 
in the framework of the one-band Eliashberg model. The effective Eliashberg function, determined 
on the basis of the transport function (R.S. Gonnelli, et a/., Physica C 341, 1957 (2000)), was used 
during calculations. It was shown that the dimensionless ratios are equal to: R± = 2 A (0) /UbTc — 
3.87, R 2 = AC (T c ) /C N (T c ) = 1.79 and R 3 = T C C N (T c ) /H 2 C (0) = 0.159. The value Rt fairly 
agrees with the experimental data whereas R2 and R3 agree very well. 
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I. INTRODUCTION 

The superconducting state in the most of the low 
temperature superconductors is being induced by the 
electron-phonon interaction [l|, Q. The coupling be- 
tween the electron and phonon system can be modeled 
with a use of the Frohlich Hamiltonian (3[. Let us no- 
tice that using a canonical transformation, during the 
elimination of the phonon degrees of freedom in Frohlich 
operator, it is possible to obtain the Hamiltonian of the 
BCS theory [J], [5j. It should be clearly marked that the 
BCS model is able to accurately describe the thermody- 
namic properties of the low temperature superconductors 
in the limit of the weak coupling, i.e. for A < 0.2, where 
A denotes the electron-phonon coupling constant. 

In order to precisely estimate the thermodynamic pa- 
rameters of the superconductors that have a greater value 
of A, the approach proposed by Eliashberg should be used 
In the Eliashberg scheme the analysis of the su- 
perconductivity issue is being started directly from the 
Frohlich Hamiltonian, which is initially written in the 
Nambu notation Q. Next, with a use of the matrix Mat- 
subara functions the Dyson equations are determined; 
the self-energy of the system is being calculated with an 
accuracy to the second order in the equations of motion 
[l|, Q, In the last step, the Eliashberg set is deter- 
mined in a self-consistent way. 

The set of the Eliashberg equations generalizes the fun- 
damental equation of the BCS model. In particular, one 
can take into consideration the complicated form of the 
electron-phonon interaction with a use of the Eliashberg 
function. Moreover, the application of the full version of 
the self-consistent method enables to estimate the elec- 
tron band effective mass in the presence of the electron- 
phonon interaction and the energy shift function that 
renormalizes the electron band energy [l|. 
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From the mathematical point of view the solution of 
the Eliashberg set is a truly complicated matter. For that 
reason, the effect of the electron band energy renormal- 
ization is usually omitted what results in the reduction 
of the equation's number by one third. The approxima- 
tion given above, in the most cases, does not affect the 
final results in a significant manner. This is evidenced 
by the good agreement between theoretical predictions 
and the experimental data Q . It need to be marked that 
the strict solution of the simplified Eliashberg set is not 
easy and can be performed only with a use of a powerful 
computer and highly advanced numerical methods [|[. 

The electron-phonon coupling constant for YM2B2C 
superconductor is equal to 0.676 [lflj. In this case, the ex- 
act estimation of the thermodynamic parameters is pos- 
sible only in the framework of the Eliashberg approach. 
In the literature related to the topic there are some re- 
ports that the simplest version of the Eliashberg equa- 
tions might not be able to cope with the description of 
all relevant physical quantities. In particular, the depen- 
dence of the upper critical field {He 2) on the temperature 
[l2| or the results obtained with a use of the direc- 
tional point-contact spectroscopy [l3j |. [lij suggest the 
necessity of the two-band model being used. On the other 
hand, some researchers basing on the thermal and spec- 
troscopic experiments are moving toward the direction of 
the one-band models with a non-trivial wave symmetry 
{s + g or even d-wave symmetry) [l5j]-[l8|. The issue is 
clouding by the fact that the calculations conducted till 
the present day for the one-band Eliashberg model were 
not strict; the study were based on very simplified form of 
the Eliashberg function [l9| or approximated analytical 
formulas 0. 

For this reson, in the presented paper, the most im- 
portant thermodynamic parameters for YM2B2C su- 
perconductor were analyzed exactly in the framework 
of the one-band Eliashberg model. The calcula- 
tions were based on the effective Eliashberg function: 
[a 2 F(^)] eff = 1.283 [a 2 F(^)] tr , where the transport 
function ([a 2 F(^)] tr ) was determined by Gonnelli et al. 
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in [19j (for details see also jlQj). 



II. THE ELIASHBERG EQUATIONS 

The Eliashberg equations on the imaginary axis can be 
written in the following form [6]: 
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where A n = A (iuj n ) represents the order parameter 
and Z n = Z (iuj n ) is the wave function renormaliza- 
tion factor; the n-th Matsubara frequency is denoted by: 
uo n = ^ (2n — 1), where j3 = l/ksT {ks is the Boltzmann 
constant). 

In the Eliashberg formalism the coupling between the 
electron and phonon system is considered with using of 
the formula: 
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The effective Eliashberg function [a 2 F (£7)] eff for 
YM2B2C superconductor is shown in Fig. [TJ On the basis 
of the presented data, it is easy to see an exceptionally 
strong coupling that appears between the electron gas 
and the crystal lattice vibrations of frequency of about 
20 and 50 meV. The mentioned effect was confirmed with 
a use of the point-contact spectroscopy measurements 
in [2l[. The value of the maximum phonon frequency 
(O max ) is equal to 67.51 meV. 

From the mathematical point of view, the strict anal- 
ysis of the Eliashberg equations on the imaginary axis, 
however complicated, is much simpler than on the real 
axis. In the imaginary domain the arguments of func- 
tions A n and Z n take discreet values, thus all problems 
connected with the numerical analysis of the non-linear 
integral equations are eliminated [6[. On the imaginary 
axis the Eliashberg equations can be solved in the follow- 
ing way [9| : in the first step one need to define functions 

= A n Z n and Z n respec- 

K (u n - uj m ) ±K (u n - cj_ m+ i). 
Next, the fact that functions (j) n and Z n are symmetric 
is used. The Eliashberg equations take the form: 
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FIG. 1: The form of the effective Eliashberg function deter- 
mined by using the transport function obtained by Gonnelli 
et al [3. 
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The symbol 5 njiri appearing in Eq. (O denotes the Kro- 
necker delta. Let us notice, that the parameter M must 
be chosen in a such way, that the solutions of the Eliash- 
berg equations for the large values of n and T = [T] min 
would take their asymptotic form: A n ~ and Z n ~ 1. 
In the case of YM2B2C it is enough to assume M = 800. 
Finally, the Eliashberg set is solved in an iterative way. 



III. THE NUMERICAL RESULTS 

The Eliashberg equations were solved for the temper- 
ature range from ks P1 m i n = 0.2 meV to ksTc = 1.335 
meV. The form of the order parameter for the selected 
values of the temperature is presented in Fig.O It is 
easy to notice that together with the decreasing of the 
temperature, the function A n takes the higher maximum 
values (always for n = 1) and it becomes wider. The de- 
pendence of the order parameter on the temperature can 
be traced, in the most convenient way, by plotting the 
function A n= i (T); see inset in Fig. [2j In the considered 
case we have taken 350 accurate numerical values of the 
order parameter for n = 1. We notice that the function 
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A n= i (T) can be fitted by the simple formula: 



A n=1 (r) = A n=1 (o)yi-(^ , (8) 

where A n= i (0) = A n=1 (k B T = 0.2 meV) = 2.559 
meV and j3 = 3.21. 
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FIG. 2: The order parameter function on the imaginary axis 
for the selected values of the temperature; the first 100 values 
of A n are presented. The dependence of the order parameter 
A n =i on the temperature is plotted in the inset. 

The familiarity with the form of the function A n for 
ksT = 0.2 meV enables the calculation of the order 
parameter value near the temperature of zero Kelvin 
(A(0)). In particular it was assumed that A (0) ~ 
A (&bT = 0.2 meV). In order to achieve that, the order 
parameter on the imaginary axis needs to be analytically 
continued on the real axis (A (a;)), deriv ing the coeffi- 
cients pAj and qAj in the expression [22| , [23j : 
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where r = 400. The plot of the real (Re) and imagi- 
nary (Im) part of the function A (uj) is shown in Fig. [3j 
In the last step, the parameter A (0) should be calcu- 
lated on the basis of the equation [2|, [6[: A (T) = 
Re [A (u = A (T) , T)]. The value 2.584 meV was ob- 
tained. 

The determination of the parameter A (0) allows to es- 
timate the dimensionless ratio: R\ = ? A ffl which, in the 

1 kbj-c ' 

BCS theory, is the universal constant of the model and 
[i?i] BCS = 3-53 [5]. In the case of YM2B2C, the greater 
value of Ri, equal to 3.87, was obtained. With reference 
to the result predicted by the two-band Eliashberg model, 
the value of R\ estimated in the framework of the one- 
band model should be interpreted as the resultant quan- 
tity. In particular the two-band Eliashberg model, which 



very well reconstructs the experimental value of R\ , pre- 
dicts: R t = 0.71[i?i] z + 0.29[i?i] s = 3.78, where [R^ 
and [Ri] s denote the ratios for the large and the small 
value of the order parameter respectively [24[. When 
comparing the obtained results it is easy to notice, that 
the one-band model insignificantly overestimates R\. 




FIG. 3: The real and imaginary part of the order parameter 
on the real axis; ksT = 0.2 meV was assumed. 

In Fig. H] the dependence of the wave function renor- 
malization factor on the Matsubara frequency for the se- 
lected values of the temperature is presented. On the 
basis of the achieved results it has been stated, that 
function Z n takes its maximum always for n = 1. In 
the framework of the Eliashberg formalism the value 
Z n= i plays very important role, because it determines 
the ratio: Z n= \ = m* e /m e , where m* denotes the elec- 
tron band effective mass in the presence of the electron- 
phonon coupling and m e is the electron band effective 
mass in absence of the electron-phonon interaction. Bas- 
ing on the results presented in Figure's |4] inset, it has 
been concluded, that m* takes its maximum value equal 
to 1.676ra e for T = Tc- Let us mark the fact, that in 
the case T = Tc, the electron band effective mass in the 
presence of the electron-phonon coupling can be calcu- 
lated with an use of the formula: 

m* = (1 + A)m e . (10) 

On the basis of expression ([TO]) it has been stated, that 
the calculated value of m* is identical with the value ob- 
tained using the Eliashberg equations. The above result 
partially confirms the accurateness of the advanced nu- 
merical calculations. 

Next, we have calculated the following ratios: 

AC(T C ) 



R2 



C N (T C ) 



and 



= T C C N (T c ) 
~ H c (0) 



(11) 



(12) 



4 



N 



1.7 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 
1.0 



k B T=0.2 meV 
k B T=0.4 meV 
k B T=0.7 meV 
k B T=1.0 meV 
k B T c =1.335 meV^ 



1.67 



M.66 



1.65 




kj meV 

D 





0.3 0.6 0.9 1.2 



k B T meV 



C7k B p(°) meV 

C N /kn(°) me W 



_i i i_ 



0.4 0.8 1.2 1.6 
k J meV 

D 



FIG. 4: The wave function renormalization factor for the se- 
lected values of the temperature; the first 400 values of Z n 
are presented. The dependence of the wave function renor- 
malization factor on the temperature for the first Matsubara 
frequency is plotted in the inset. 



FIG. 5: (A) The dependence of the free energy difference 
between the superconducting and normal state on the tem- 
perature. (B) The specific heat of the superconducting and 
normal state as a function of the temperature. 



where AC = C s — C N stands for the difference between 
specific heat of the superconducting and normal state. 
The dependence of AC on the temperature is being de- 
termined with a use of an expression: 



AC (T) _^ 1 d 2 [AF/p ( 
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In Eq. (fT3|) the symbol AF denotes the difference be- 
tween the free energy of the superconducting and normal 
state; p (0) is the value of the electron density of states 
at the Fermi level. 

On the other hand, the specific heat in the normal state 



can be calculated using the formula: 
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field (He) should be estimated in accordance with the 
expression: Hc — 



Having an open solutions of the Eliashberg equations, 
AF is determined directly from |25j : 
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where symbols and denote the wave function 
renormalization factor for the superconducting state and 
the normal state respectively. 

The determined form of AF (T) is presented in Fig. 
[5] (A), while in Fig. 0(B) we are shown the dependence 
of the specific heat on the temperature for the supercon- 
ducting and normal state. The characteristic jump of the 



specific heat that appears at the critical temperature was 
marked with a vertical line. We notice that the specific 
heat of the superconducting state was obtained on the 
basis of the 350 accurate values of AF/p (0). 

On the basis of the determined thermodynamic func- 
tions, the parameters R2 and Rs were calculated. In the 
case of YM2B2C these coefficients take following values: 
1.79 and 0.159. Let us notice, that R2 and Rs are the uni- 
versal constants of the BCS model and [^2]bcs = 1-43 
and [Rs] BCS = 0.168 0. When basing on the above 
results one can clearly see, that for YM2B2C supercon- 
ductor the value of ratio R2 significantly deviates from 
the BCS prediction. 

The experimental ratios R2 and Rs were determined in 
[20| . The following results were obtained: [^2] exp = 1.77 
and [i?3] exp = 0.160. When comparing our theoretical 
values of R2 and Rs with the experimental ones, it can be 
easily noticed that the one-band Eliashberg model prop- 
erly describes the experimental data. 



IV. CONCLUDING REMARKS 

In the framework of the one-band Eliashberg model 
the selected thermodynamic properties of YM2B2C su- 
perconductor was calculated. In particular, the funda- 
mental ratios R\ ~ Rs were determined. It has been 
stated that, in the case of R\ the one-band Eliashberg 
model predicts a value non-significantly higher than the 
value determined by the two-band model, which deter- 
mines R\ with a very high accuracy. The result above 
is connected with the simplified description of the super- 
conducting phase in the one-band model in comparison 
with the two-band model; we have taken the one effec- 
tive Eliashberg function instead of the four Eliashberg 
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functions and the four elements of the Coulomb pseu- 
dopotential matrix. 

In the cases of the two remain ratios - the agreement 
between predictions of the one-band Eliashberg model 
and the experimental data is very good. It is worth to 
notice that, in the opposition to R\ and i?3, the value of 
R2 significantly deviates from the prediction of the BCS 
model. 

In the last part of the summary let us turn a reader's 
attention toward the problem of the correct determina- 
tion of the upper critical field. When basing on the 
papers cited in the introduction one can suppose with 
a large probability, that the exact form of the function 
Hc2 (T) can be reproduced only in the framework of the 
two-band theory. Let us remind that the two-band ap- 
proach was successfully used by us for description of the 
thermodynamic properties of the superconducting state 



inducing in MgB2 [26]. However, in the case of YM2B2C 
discussed problem is far more complicated because the 
appropriate Eliashberg functions have not been calcu- 
lated in the branch press (only electron-phonon coupling 
constants are being used in the existing two-band ap- 
proach) . At present this issue is being intensively studied 
by us with a use of the ab initio approach [271 ]. 
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